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Reduced Model Number cf States 

White Process Noise Ccvariance Matrix (Ng,Ng) 
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Measurement Vector (Nc, 1) 


ACKNOWLEDGEMENTS 


I would like to express my sincere appreciation to 
Professor D.J. Collins whose assistance and encouragement 
enabled this research. I would also like to thank Professor 
D.G. Williams, Director of the W.R. Church Computer Center, 
for his assistance and cooperation in providing the addi- 
tional computer memory allocation which was essential to the 
completion of this work. 

This thesis is dedicated to my wife Gayle, my daughter 
Kristin and my Son » ‘Scott: Without their constant love, 
support and understanding this work would not have been 
possibtle. 


10 


I. INTRODUCTION 

Modern control methods allow the control engineer to 
develop a control system to regulate every parameter ofa 
physical systen. He has been aided by computer simulations 
of non-linear systems and linear approximations of those 
systems. The desired contrcl system can be developed 
parallel to or even ahead of the physical system which is to 
be controlled. 

In ,many design applications areas, the engineer is 
confronted with designing a ccntrol system which can be 
represented analytically by a very large number of first 
order, linear, time invariant differential equations. This 
is especially true with highly complex and aero-elastic 
aircraft and space vehicles. They analytical “model Ws” a 
prime ingredient in the controller design process for any 
design technigue. 

Many practical limitations exist, such as computational 
requirements, that require the order (number of equations) 
of this complex model to be reduced for analysis, synthesis 
and actual implementation of the control system. Some method 
is required to analyze these dynamical models and establish 
Simpler systems which include elements critical to the 
desired control function. Without’ ’such “sSzmplifi'cati on, 
application’ of design “procedures” can result “in highly 
complex and parameter-sensitive controlled systems. With a 
reduced order model, regulator synthesis procedures become 
more intuitive and far less sensitive to parameter varia- 
tvon . The reduced order model is far easier to handle in 
formulating the control system but must characterize the 
physical system with sufficient accuracy such that perforn- 
ance objectives for the controlled physical system can be 


met by designing control laws from the reduced order model. 
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The intent of this thesis is to present the mathematical 
basis for creating such a reduced order model and to develop 
the actual computer software which creates that reduced 
model based on a given large order system. This computer 
program (OPTRED) is developed to interface with an existing 
control analysis FORTRAN program named OPTSYS. This thesis 
will discuss the modifications or additions to the iPUS¥e 
program which are necessary to that interface. 

The OPTimal SYStems control . program was oOriginaley 
developed by Hall [Ref. 1] for the study and application o£ 
optimal systems control theory. Later modifications tc the 
program were made by Walker [Ref. 2] and liu [Ref= 3i§ae 
Stanford University. Hoden [Ref. 4] modified the OPTSYS 
program to present a user-friendly, interactive version of 
the. progzran, , (OPTSYSx)\- Diel [Ref. 5] introduced changes 
which enabled the user to save various data for re-entry 
into the program and also to create data files which were 
necessary for the execution of atime response program 
{OPTCALC). Further modificaticns have been made by Laptas 
[Ref. 6] which were necessary tc create data sets for input 
to his OPTGRAPH program. This FORTRAN program enabled the 
user to obtain Pole-Zero, Root-Locus, Nyquist, Bode and 
Nichols. plots. 

It is assumed that the user has a basic understanding of 
the fundamental concepts of control theory and optimal 
systems design. The symbology conventions... of. BEyson 
[Ref. 7] are used in the discussion of program operation and 
system descriptions. 

An overview of the OPTSYSX progran, its capabilities, 
and the modifications to that program which are necessary 
for the interface with OPTRED is presented first. This is 
followed by the mathematical tasis upon which OPTRED is 
formulated and by a full description of the operation of 


that progran. Finally, examples of large order control 


2 


eee a cA RR a 


systems and the corresponding reduced models which OPTRED 
generates are presented and analyzed. Those systems include 
a sixteenth order model of the .F100 turbofan engine anda 
ninety-eighth order model of the X-29A aircraft. Complete 
PeegGbam Jastangs for “the OPTSYSX “program and the OPTRED 


program are included as Appendices A and B respectively. 
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II. THE OPTSYSX COMPUTER PROGRAM 


A. PROGRAM OVERVIEW 


OFTSYSX «is “ajdoublesy precision; interactive FORTRAN 
program which employs modern ccntrol theory analysis tech- 
nigues. It is developed to be compiled and executed by the 
Naval Postgraduate School's IBM 3033 System 360/370. tes 
primary capabilities include the calculation of the open 
loop eigensysten, and the fixed closed loop system; the 
synthesis of optimal regulators or filters; the power 
spectral density, and modal distribution computations. 

The fundamental system equations used by the OPTSYSX 
program for its computations are of the state variable forn. 


The basic system equations are: 


system model 

x = [F]*x + [G]*u + [GAM]*w (2a 
output equation 

y = LH jtx ss LoaFu (2.2) 
measurement equaticn 

z= [H]*x + v (2-3) 
estimator equation 


xe = [F]*xe + [G]*u + (K]*(z - [H]*xe) (2.4) 
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open loop transfer function 

PaaS Me) Sere |=241 Gy (2.5) 
closed loop noise transfer function 

[H}*{/s[I] - [FJ]]-'*[ Gan] (2.56) 


compensator transfer function from measurement to input 


Peg Siri Gtr ed si Ke P ap et+Lk | (2.7) 
where 

x = state vector (Ns X 1) 

x = derivative of the state vector (Ns X 1) 

u =“control*vector (Ne “xX ‘1) 

y = output vector (No X 1) 

Zz, = measurement vector (No X 1) 

W = white process noise vector (Ng X 1) 

[Pra = open-loop dynamics or plant matrix (Ns X Ns) 
{G] = control distributicn matrix (Ns X Nc) 


[Gam] = state disturbance distribution matrix (Ns X Ng) 


[H] = measurement distribution matrix (No X Ns) 
[D4 = control feed-forward matrix (No X Nc) 

Pes = control gain matrix (Nc X Ns) 

bE = identity matrix (Ns X Ns) 


OPTSYSX is an extremely large and complex program which 
contains over 4000 lines of code. LES (existing jstandard 
dimensions will accomodate a thirty-second order system with 
up to 10 controls and 10 observations or measurements. TO 
execute even this version of OFTSYSX the user requires one 


megabyte of virtual storage. Any “Lurther | signafiicant 


15 


additions to that program's size would cause user difficulty 
in obtaining sufficient memory for program execution. It is 
for this reason that Diel"s OPTICALC, Laptas"*. OPTCGRAPH Gane 
the program developed by this thesis have been developed as 
independent programs which depend upon various data created 
by ODTSYSK. The minor modifications this thesis introduces 
to the OPTSYSX program are not detrimental to its original 
capabilities. 


Be. OPTSYSX MODIFICATIONS 


This thesis is ccnfronted with the problem of analyzing 
very large control system models. Therefore, the additions 
of Hoden's subroutine SETUP and Diel's subroutines RDMATF, 
RDMAT, and WRITMAT were crucial to system analysis. These 
routines make multiple computer analysis runs’ possible 
without painstaking re-entry of system parameters and 
matrices. In particular,,; the system "E";,-"C", oes 
"Gam" are some or the matrices which are stored for re-entry 
into OPTSYSK or stored as data for use by, otheryanalyems 
programs. Subroutine WRIMAT is responsible for the storage 
of these matrices inthe data file "OPTMAT". Its jis thas 
file which OPTRED utilizes. te input all pertinent Sea 
system data. 

In many control applications, a control feed-forward 
matrix ("D") is present in the system model. Therefore, 
this thesis makes provisions within OPTSYSX which enable the 
user to save the "D" matrix in the same manner by which the 
other system matrices are _ saved. Subroutine SETUP, whaer 
enables the user to input system matrices from file data, 
was changed to enable the user to input the "D" matrix. 
This consisted of adding "D" to the routine's parameter list 
and also added a read statement to input that matrix froma 


data file. Present routines in OPYSYSX which "save" the 
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oe 


wpm ong, 6AM, 6©6h6Uand "GAM" matrices were complemented by the 
addition of routines which also save the "D" matrix. These 
hew routines are analagous to existing "save" routines. In 
addition to these routines a new flag, ISAD, was added to 
miermexasting flags TSAM, \/fSAGC,s ISAH, “and IGAM. This new 
flag provides the same function in "saving" matrices as do 
the existing flags. For program continuity, it was also 
mecessary to add ISAD to ‘subreutine INNER, which is the 
driving subroutine in the OPTSYSX progran. 

Subroutines RDMATF, RDMAT, and WRITMAT also reguired 
minor changes to accomodate the "D" matrix. These changes 
consisted of the addition of "D" to the subroutine parameter 
lists and the addition of read or write statements for 
Palermo c-onthy- sor "Dante ,OPISYSX or for writing ‘’-"D" to 
the OPTMAT data file. The exastang flag, LFDFW, was also 
added to the parameter lists of these routines. Wheres: pliag 
indicates whether or not a "D" matrix will be input by the 
user. OPTRED uses this flag in its computational methods in 
the same manner as does OPTSYSX. Appendix A lists’ the 
OPTSYSX program as modified by this thesis. 
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IIIl. ANALYTICAL DEVELOPMENT OF A REDUCED ORDER MODEL 

There are numerous methods of order reduction for 
control systems as discussed by Enns [Ref. 8] of Stanford 
University. This thesis uses the theory discussed by DeHoff 
and Hall [Ref. 9] which presents a method of reducing the 
order of the linear system based on dominant states chosen 
by the designer. This dominance is determined by modal 
analysis of the eigensystem under investigation. Dougherty 
[Ref. 10] utilized this same methodology in the development 
of a computational model used in the analysis of the F100 
turbofan engine. A summary of this technigue is provided 
below. 

The method of reduction involves finding the eigenvector 
transforn matrix, T,. .such. that 


X*¥T = TA (32%) 


where A is the diagonal eigenvalue matrix. Then, by 


defining an alternate state vector, Z, where 


X = Tz (3.2) 


the original system may be re-written in modal coordinates 


as 


Z = A¥Z + T-1*G*U (3.3) 
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The eigenvalues which are tc be included in the initial 
reduced order model are determined through modal analysis. 
The eigenvalue matrix can be tre-ordered by constructing a 
Matrix Q, which has the value 1 in the position corre- 
sponding to the eigenvalue to te included (column) and the 
re-ordered position (row) of that eigenvalue. For example, 
ma cigenvalue . 20). of the original, A.) matrix, .was.,.to,..be 
re-ordered to the first position, a 1 would be placed in the 


(1,20) position of the Q matrix. Now, the reordered matrix 


At = C¥A (3.4) 


is incorporated in the linear model 


O¥Z = A'KZ + Q¥T-1%G4U (3.5) 
letting 
Z" = Q*Z (3.6) 
thus 
ZY = A'*Q-147" + Q¥T-14G%U (32) 
and 
X = T#Q-1%z' (3.8) 


A Matrix R can be constructed such that 


NS) 


fe 
X' = R¥X = T¥Z' = |---- (3.9) 
xX 
2 
This “X' ‘vector is “"“the reordered "state weetor, where the 


elements contained in the subvector Xi are those states 
associated with the designers choice of eigenvalues while 
the subvector xX, contains all other states. 

The linear system may now be rewritten 


KX" = R¥PeR-1*X"' + R¥GHT (3.10) 
and 
Y = H¥R-1*X" + D¥*U (3. ta) 
Now define 
F' = R*¥F4R-1 (3.12) 
G' = R*G (3-13) 
H' = H*R-1 (3.14) 
and 
F F 
thay 3p Ly 
Ft = rai (32 15} 
| 
eee ee 


G? = —— 


H' = [H 


If the states in Xy tpuly.nodel the. full order 


then the subvector X, will be essentially zero and 


2 


This is incorporated into Equation 3.18 and 


re ee ee 2 As ae TWIN (CR a Spike 3 2 ame 5 Ce) 
a ( Tle) 422 21) al ( eye a2) 22 2) 


and the output equation becomes 


io= (hen = -1% re HY of ae 
oe opm qe ge aha 5 eo) 


The reduced order model can now be written as 


X = Forex. + Sr*u 
1 i 
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(3.16) 


(3547) 


systen, 


(3.18) 


(3.19) 


(8220) 


(22 


(3522) 


(37523) 


Y CSP HEF. ot 
i 


where 


Gr = Gy — e/a *F , 


BES =o 
1 


DE = 7D) = VEE 


Dr *U (3.24) 


Fo, 2m aa om (sa 
1) Soo ne ee 
HS*P Gott (3.27) 
SES eee (3.28) 


This is the mathematical prccedure which OPTRED uses in 


computing reduced order models. 


is included as Appendix B. 
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The complete OPTRED listing 


Iv. THE OPTRED PROGRAM 


A. PROGRAM OVERVIEW 


OPTRED (iS an interactive, double precision FORTRAN 
mprogram which uses.a specifically formatted data file to 
apput full system data. The required matrix inputs fron 
ees, file are the full system. "F", GUE He anid pt 
Matrices. The system flags and parameters which OPTRED must 
also read for program execution are Ns, Nc, No, and IFDFW. 
The OPTSYSX program creates OPTMAT data which is the source 
of this required data. Through interactive input, the user 
specifies all other data which is required to create reduced 
model data. 

OPTRED also reads from the CPTMAT file the noise, feed- 
Pack) control and-~“the “output and control cost weighting 
matrices. These are not required inputs but only input to 
the program to facilitate creation of a reduced order data 
file which is identical in format to the OPTMAT data file. 
The final result of OPTRED is the creation of OPTMATR, a 
Gata’ file which contains the reduced system "F", “Gu, uy 
ena o"D"ematrices’ and iis” formatted for re-entry into the 
OPTSYSX progran. 

heEaS possible for the usier to construct*a“VYdata file 
hamed OPTMAT to provide input data to OPTRED. However, only 
in rare applications would the full system not be analyzed 
prior to the development of the reduced nodel. OPTSY¥SK is 
therefore the simplest and most preferred method of data 
generation. 

ZO) execute the  OPTRED prcgran, a minimum of three 


separate files must be immediately accessible to the user: 
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1) The compiled text file of OPTRED 
2) The OPTMAT data file 
3) The OPTRED EXEC file 
The first two required files are self explanatory. The 
OPTRED EXEC file sets up the required read/write file defi- 
nitions and calls the necessary library functions for execu- 
tion using the IBM 3033 Systen. This EXEC allows the user 
to execute the OPTRED program by simply entering the werd 
“NOPTRED" at the terminal. The program then continues iD a 
user-interactive mode until completion as described in the 
next section. The OPTRED EXEC file is listed in Appendage 


B. PROGRAM OPERATION 


OPTRED was written to satisfy two specific objectives. 
The first is to provide a methcd to assist in the analysis 
of actual. lagge: order,;control. isystenss The seconds» ise 
create a program to be used aS an instructional .toolger 
students studying controls related problems. 

To accomplish these objectives the emphasis while 
creating OPTRED was to make the program as user oriented as 
possible by eliminating data entry ambiguities and providing 
features which prevent the user from other inadvertant or 
invalid program entries. A serious effort was also made to 
Minimize the amount of virtual memory required to execute 


this frogran. 
and Flow 


OPTRED is composed of one main program and _ seven 
subroutines. The program can be divided into three kasic 
catagories: 

1) "Bale Data,Input 
2), interactive Data .nput 
3) Calculation 
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Both file and interactive data entry is accomplished 
by the main program and by the primary subroutine, REDUCX. 
She) actual computation Yor the reduced order” model ‘is 
performed by subroutine REDUCX which utilizes subroutines MI 
amas i MAMULT for itsimatrix  lnversion sand ,multiplication 
requirements. The remaining subroutines) jare “utziized 
throughout OPTRED for interactive entry of numerical -and 
character string data. The program requires no external 
subroutines for its execution. 

The main program begins by presenting the user with 
a brief operational description of OPTRED. At this point, 
the user is given the option of having the program's general 
results repeated to the screen cr printed to a listing file. 
OPTRED creates a permanent data file which contains the 
reduced data regardless of the user's option here. The 
program then reads full system, flags and “parameters, 
excluding matrices, from the OPTMAT data file. A brief 
description of the full system under investigation is then 
presented as an indication to the user that the desired data 
as ebeing @nput. This description includes the order of the 
system, the number of controls and observations and whether 
Seenot aysysten "D" matrix, willbe input. The latter is 
indicated either by " Wor "net" in the following example < 
The program then prcmpts the user for entry of the desired 
order of the reduced systen. 

The main prcegram presentation as it appears during 
execution: 

OPTRED 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
SEG" MH", AND "D"™ (IF INPUT) MATRICES. 


THE FULL SYSTEM MATRICES MUST BE READ FROM A FILE NAMED 
"OPTMAT DATA" AS CREATED BY THE OPTSYSX PROGRAM. 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
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AND THE ACTUAL STATE #"S (IN ASCENDING ORDER) WHICH REPRESENT 
THE REDUCED MODEL. THE ORDER OF THE REDUCED ORDER MODEL MUST 
BE LESS THAN THE ORDER OF THE FULL SYSTEM. 


DO YOU WISH TO CONTINUE? 
LY PE VYES™ ORR NON. 
DO YOU WISH RESULTS TO SCREEN O8 DISK? 


NOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARILESS OF YOUR ANSWER. 


TYPE "S™ FOR’ SCREEN OR "“D(FOR Drsk. 


THE ORDER OF THE FULL SYSTEM IS; 
THE NUMBER OF CONTROLS IS: 

THE NUMBER OF OBSERVATIONS IS: 

A "D" MATRIX WILL BE INPUT. 


DO YOU WISH TO CONTINUE? 
TYPES VYES “OR ROW. 


ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. 


Subroutine REDUCX is now called and immediately 
reads the full system matrices from file data. The user is 
then prompted for entry of the significant states which will 
compose his reduced model. This program makes a provision 
for entry of significant states? £eom a specifically 
formatted, user created data file. At this point in gene 
program he is given that option. Once the significant 
states have been entered, the program repeats those states 
to the screen and the user is given the option to make 
changes at that time. If the significant states are satis- 
factory to the user, the program continues with the computa- 
tion of the reduced order model. At the completion "ez 
OPTRED, the reduced system data is created and can be exan- 
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med tfochr stably, controllakvlity “and observability by 
Gener andlysis techniques aS discussed in the following 
section. Prior t@, exiting the OPTRED progran, the user is 
peovided the option of re-executing OPTRED if he finds 
mtttal results unsatisfactory. A typical recording session 
mech deprcts this program filcwis presented in the next 


SeCt iON. 
Z2syeUser ‘Protection Features 


The user of OPTRED is protected from abnormal 
program termination in several ways. Initially, the program 
presents its capabilities and user required entries. If he 
aS not prepared to continue at this point, the user is 
offered the option of exiting the progran. Next,| the fuadl 
system parameters are presented and the user is once again 
able to exit if these parameters are unexpected. 

OPTRED requires that the order of the reduced systen 
bewgreater than zerc but less than the order of the full 


systen. If the user inadvertantly enters a number which is 
our Of this range, the program issues a warning to that 
effect and allows recovery. The computation of the reduced 


order model relies ufon the construction of a reduced state 
Matrix. OPTRED requires that these states be entered in the 
order that they would normally appear in this matrix. For 
fats reason, the user must enter the significant state 
Numbers in ascending crder and those state numbers must lie 
in the range from one to the frumber of states in the full 
systen. If the user inadvertantly enters these state 
humbers out of order or out of range the program frompts a 
warning and allows recovery. 

The method of reduced order computation involves the 
inversion of the state submatrix P50! as described in 
Siapter Til, Equations’3.25 threugh 3.28.' The nature of the 


plant matrix and the selection of reduced order states can 
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be coupled such that this matrix will be singular ~or Mens 
invertable. OPTRED will detect this condition and issue 
information to the user that a reduced order model cannot be 
calculated for that full system from the given desived 
reduced states. Reduced data will still be computed but the 
user is warned that this data is invalid. 

Subroutines RDINT, RDREAL and RDCHAR are responsible 
for the interactive input of data and expect an integer 
input, a real number input and a‘ logical "YES" o&e"NO™ Spee 
respectively. If the user inadvertantly makes an incompat- 
ible entry these subroutines issue warnings and allow 
another opportunity for data entry. The entry of a null 
line is included in these improper actions and the entry of 
two consecutive null strings will cause termination of the 
progran. This function allows the user a further means of 


exiting OPTRED if he so desires. 
3. Large Order Systems 


The analysis of large crder systems, particularly 
the ninety-eigth order model, presents a major problem to 
the user with limited virtual memory assets. The longitu- 
dinal control system analyzed by this thesis has a (98 X 98) 
"et matrix, a (2 X 98) "H" matrix and a (98 X 1) “"G" nates 
To analyze this systen, OPTSYSX must be executed in its 
"increased dimension" form. Io compile and execute this 
version of OPTSYSX requires 2.5 Megabytes of virtual memory. 

Both OPTRED and OPTGRAPH are also dimensioned to 
enable analysis of a ninety-eigth order model. When 
OPTSYSX, OPTRED, OPTGRAPH and their peripheral data are all 
used during a typical analysis session, very careful atten- 
tion to file and virtual memory management is mandatory to 
prevent the user from exceeding his virtual machine's 


capacity. 
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C. INTERFACE WITH OTHER ANALYSIS PROGRAMS 


To sufficiently analyze a control system it is desirable 
EOwucrlizelOPTSYSX,;, “OPTCALC, CPTGRAPH and OPTRED. It may 
also be desirable to conduct this analysis during a single, 
continuous interactive computer Session. However, in the 
analysis of very large systems this concept presents serious 
difficulties in the areas of program flow, data management 
and physical memory assets. For this reason, OPTRED has not 
been developed to be fully autcmatic in its interface with 
these other analysis programs. 

The OPTRED program reguires the existance of the OPTMAT 
data file. As previously discussed, OPTRED generates the 
OPTMATR data file which contains the reduced model data. To 
analyze this data using OPTSYSX, the user must now rename 
the OFTMATR data file for entry into OPTSYSX as OPTMAT data. 
Care must be exercised at this point to either erase the 
full system data file or to rename it to prevent ambiguous 
data files. Although this procedure may seen cumbersome at 
first, it is easily accomplished in the XEDIT mode and 
provides the user with the capakility of maintaining several 
system data files which are readily available for analysis. 

After OPTSYSX processes the reduced data, additional 
data files are created which enable the execution of other 
system analysis programs. OPTGRAPH and OPTCALC utilize the 
OPTMAT and OPTGROL data files, respectively, to provide the 
functions described in chapter 2. When the user is investi- 
gating several systems it is both prudent and convenient to 
rename these files after their use for the purpose of future 


analysis. 


D. EXAMPLES OF ORDER REDUCTION 


Two large order systems will be presented to evaluate 


reduced models generated by OPTRED. ThelELTESt ‘system is a 


2:9 


sixteenth order model of the F100 turbofan engine. The 
final example is a ninety- eighth order model of the X-29A 
aircraft's longitudinal control systea. Bode analysis was 
conducted for the full and reduced models and graphical 
comparisons follow at the end of the chapter. Recorded 
terminal sessions for each example will be presented to 


fully illustrate actual program operation. 


1. Example of a Sixteenth Order System 


The F100 turbofan engine was chosen for ainderad 
analysis and [Ref. 9] describes this system in detail. The 
systen “hwy "Gr, HY and "D" matrices “were obtained ren 
[Ref. 9: pp. 83-85] and the method of data entry to 02957 
is depicted in subroutine SETUF found in Appendix A. The 
selection of significant states is based upon a desired 
control bandwidth of 1 to 10 Hertz and following modal anal- 
ysis these states were chosen as 1, 2, 5, 1l°and” tez The 
following is a computer terminal system in which a’ feeeg 


order model is generated. 


BEGIN RECORDING OF TERMINAL SESSION 
Rs .T=0.0170.03, 20305288 
OPTRED 
FILEDEF 05 TERM 
FILEDEF 03 DISK REDUCI DATA Al 
FILEDEF 06 DISK OPTRED LISTING Al 
FILEDEF (07 DISK«STA TES, DATA.A’ 
FILEDEF 09 DISK OPTMAT DATA Al 
FILEDEF 10 DISK OPTMATR DATA Al 
GLOBAL TXTLIB FORTMOD2Z MOD2EEH IMSLDP NONIMSL 
LOAD OPTRED ( START 
EXECUTION BEGINS... 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
tee wen ve, AND "DY (IF INPUT) ARR ECE. 
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Te ULL SvoLe MN MATRTCESUMUST ABE READ FROM, .A. FILE NAMED 
POPTMAT DATAY AS CREATED(BY THE OPTSYSX PROGRAM. 


yOU MUST ALSO, ENTER THE DESIRED REDUCED ORDER (NUMBER OF. STATES) 
arp THE ACTUAL STATE #"S (IN ASCENDING ORDER) WHICH REPRESENT 
Tee REDUCED MODEL. THE ORDER OF THE REDUCED MODEL MUST BE 

feos THAN THE ORDER OF THE FULL SYSTEM. 


DO YOU WELSH TO CONTINUE? 


EYP RULES OR UNO. 


YES 
DO YOU WISH RESULTS, TO SCKEEN OR DISK? 


NOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARDLESS OF YOUR ANSWER. 


TYRE UST FORGSCREENWCR "DY" FOR DESK. 


) 
ZHEVORDER OF THE FULL |SYSTEM,» ES; . 16 
ZRSONUMBERIOFACCNEROLS TSs 55 
THE NUMBER OF OBSERVATIONS IS: | 
AOD MARE XW BE INPUT. 
DO YOU WESH TO CONTINUE? 
BYES UVES MOR Git. 
YES 
PNTER THE DESTREDTREDUCED VORDER OF \THE, "FF" MATRIX. 
% 
18 


**e%x%* WARNING: REDUCED ORDER MUST BE GREATER THAN 0***** 
AND LESS THAN 16 


ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. 
% 
5 
DO-YOU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED ORDER 
MODEL FROM A DATA FILE? 
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DATA FILE MUST BE NAMED "STATES DATA AI" IN FIXED SFORROae 


THE READ FORMAT SE See ioe 


PER 72 CHARACTER LINES 


SYP EUV ESY COR MON Gers 


NO 


ENTER THE "N" SIGNIFICANT STATES WHICH REPRESENT THE RElGeee 


MODEL. ENTER STATE #°S IN ASCENLTING ORDER: 


ARES 


DO YOU WISH TO CHANGE ANY OF THE SIGNIFICANT STAWESe 


TO BE CHANGED. 


STATE ete jlg= 
a 
1 
SLALe i 2 = 
Es 
2 
STATE 403 = 
fe 
3 
STATE # 4 = 
S 
11 
STATE. 4). 2= 
Re 
16 
THE REDUCED MODEL STATES 
1 2 3 11 16 
TYRE “YES ORM UNO. 
YES 
ENTER THE N= TH STATE. 
ie 
3 
ENTER NEWS SPATE #53 
5 
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THE REDUCED MODEL STATES ARE: 


1 2 


5 iby. 


16 


DO YOU GESH TO CHANGE ANY OF THE STGNEFICANT STATES? 


NO 


TP WES IUPOR LINO. 


(OUTPUT TO THE SCREEN FOLLOWS) 


THE DESIRED!) REDUCED ORDER, TS: 5 


THE REDUCED MODEL STATES ARE: 


1 2 


THE REDUCED 


=3—0 509D+00 
1.2243D-01 
Zt 2 O6D—03 
4.1314D+00 
-1.0770D+00 


THE RELUCED 


—-3.5445D-02 
-4.4094D-01 
ie 7 >26D—-02 
1.9802D+01 
—1.8909D+00 


THE REDUCED 


292381 D-01 
6.4607D-03 
0.0 
fe92535D-05 
=4.8870D-05 
fo 598D-05 
7.2415D-06 


5 1 


PLANT MATRIX 


2-4726D+00 
= P5627 1D+00 
=20 8339,1D-—03 
—5-3894)+00 
2a 2320000 


CONTROL DISTRIBUTION MATRIX 


-1.4080D+02 
=e f220DtOI 
-3.7646D+01 
3s 1639 D+02 
7.2700D+01 
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("FN") IS: 


3 9's SED OZ 
4.1060D+01 
-7.4276D+00 
= 1002.07 D*03 
1.146€7D+03 


=e 33 CIDt0 1 
7.83 06D+00 
To2 t3sp-—07 
4.9385D+01 

—3'. OF S6DF 00 


8.6913D+00 
4.2898D+00 
= te 1695D—02 
soe SODSD +0 1 
1s 9087 D£01 


("G") IS: 


Zeso 97D tON 
-1.0960D+01 
Zils S4O5D—02 
=e OO1SDtT 01 

2. 8867D+01 


OUTPUT, DISTRIBULION AATREX CH") TS 


1.0750D-01 
-9.9442D-07 
0.0 
1'3502)-05 
1.4796 D-04 
=I 0S 59D—05 
Om Zoow— 07, 


87 53D+40 2 
Oe ous D035 

0.0 
Serie Cop O02 

1.01€0D-03 
Sie) 163 D—0'2 
=i O 493 D—'02 
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1.2088D+90 
ye O22) 0S 

1.0000D+00 
"rer 39 00D 04 
De LOGlD—05 
“7 06D 04 
-—4. 9442D-05 


Jewle32D—02 
37.0229)D—0'1 
-9.4624D- 03 
= eho Dt 0) 
-4.8697D+01 


—1.8379D+04 
=o al iD +03 
mle to SLD O02 
4.8330D+04 
¢. 0035) #03 


7.0796D-02 
7 oO 05 

0.0 
-4.8476D-05 

1.4491D-05 
Poe Io OD U5 
=. 6607D—05 


THE REDUCED FEEDFORWARD MATRIX 


ey OS DOM 
=. 09 290-04 
0.0 
8. 6209D-05 
mile D0 718D--05 
1.0794D-04 
26.9257 D=05 


3.4001D+00 
=26 2600S O0n 
0.0 
=8553 10D-03 
4. 19320-0383 
5.163 96D=03 
9.8556 D-04 


(1p) 


1.0378D+01 
3.1109D-01 
0.0 
-4.4940D-03 
~9-3797D—04 
-4.0947D-04 
3-27 710-08 


IESS 


1.7019D+00 
4.8490D-03 
0.0 
GSS 92D—05 
=3'. 9608-05 
—7. I96738D—08 
—9. 4642D-05 


-4.1189D+03 
-1. 3105D+#07 
0.0 
-4. 45 130-09 
9.68:95D-—07 
-3.1434D-01 


-8.3787D-02@ 


ANALYSIS COMPLETE... YOUR REDUCED SYSTEM DATA HAS BEEN SAVED 


NO 


EPS GES 5 


IN A FILE NAMED 


NO PTMATR pDA DA. 


DO YOU DESIRE ANOTHER RUN? 


THE RESULTS OF THIS RUN WILL 


REPLACE YOUR LAST OPTMATR DATA. 


TYPE “YES 


Re T=0..3870.77 21206256 


RECORD OFF 


END RECORDING OF TERMINAL SESSION 


CR WNON. 


This is a typical computer terminal session in which 


the user requested output to 


some of the user protection 


during this session. 


Following this analysis, 
OPTMAT data 


reduced model data. 


renamed 


able characteristics, including stabubuny, 


to the full system eigenvalues. 


factory 


and further 


and 


the 
OPTESVS'X 


comparison of 


34 


the screen. 


Note also, 


that 


features of OPTRED are involked 


OPTMATR data file was 


was executed 


using 


and 


the 


The eigensystem was analyzed for desir- 
and was compared 
This analysis proved satis-— 
the: iu 


reduced 


models was accomplished by obtaining Bode diagrams for the 
system open loop transfer functions. Figures, 4.1). and 4.2 
illustrate the full system and reduced system Bode plots 
respectively. The similarity, especially within the control 
bandwidth, is readily apparent. The frequency responses 
shown are for system input 1 and output 1. Analysis of 
other input/output combinations revealed the same favorable 


comparisons. 


2. Example of a Ninety-eighth Order System 


The X-29A control system includes a ninety-eighth 
order dynamics matrix. The measurement scaling matrix is 
constructed such that there are two measured outputs. These 
outputs represent the system with and without a notch 
filter, respectively. Full system matrix data as well as 
eigensysten data was obtained for this model from 
NASA-Edwards. Full system idata was first entered into 
OPTSYSX using subroutine SETUP and the overall system was 
analyzed. Sixty Significant states were identified and are 
Shown in the following recorded terminal session. hess 
example also illustrates program execution where no “"p" 
Matrix is input and the reduced system states are entered by 
means of file data. During this session the user also 
chooses to have the general reduced model data printed to a 


listing file. 


BEGIN RECORDING OF TERMINAL SESSION 
fet —O0201/0.02 12523 356 
OPTRED 
PELEDEF 05 TERM 
FILEDEF 06 DISK OPTRED LISTING Al 
PELEDEF O07 DISK STATES DATA A‘ 
PELEDEF O9 DISK’ OPIMAT DATA A1 
FILEDEF 10 DISK OPTMATR DATA Al 
GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NONIMSL 


a5 


LOAD, OP BRED (5 START 
EXECUTION BEGINS... 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
tee Wet CA, AND “DM (Clr TNE US) Maher Gro 


THE FULL SYSTEM: MATRICES MUST BE READ FROM A SP2TLE NAR 
"OPTMAT DATA" AS’ CREATED, BY THE OPTSYSK, PROGRAM. 


SS — 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
AND THE ACTUAL STATE #"'S {IN ASCENDING: ORDER) WHICH REPRESS 
THE RELTUCED MODEL. THE ORDER OF THE REDUCED MODEL MUST BE 

LESS THAN THE ORDER OF THE FULL SYSTEM- 


Sh ee 


DO YOU WISH TO CCNTINUE? { 


TYPE SYES@ <ORy WN 


YES 
DO YOU WISH RESULTS TO SCREEN OR DISK? 
NCTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARDLESS OF YOUR ANSWER. | 
TYPE "S" FOR SCREEN OR "D" FOR DISK. 
D 
THE ORDER OF THE FULL SYSTEM IS: 98 
THE NUMBER OF CONTROLS IS: 1 
THE NUMBER OF OBSERVATIONS IS: 2 
A "D" MATRIX WILL NOT BE INPUT. 
DO YOU WISH TO CONTINUE? 
TYPE "YES" OR "NO". 
YES | 
ENTER THE DESIRED REDUCED CRDER OF THE "F" MATRIX. | 
60 


DO YOU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED ORDER 
MODEL FROM A DATA FILE? 
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DATA EEL, MUST, BE NAMED “STATES DATA. AU IN FIXED, FORMA 
Hae RAD FORMAT GS “iSr5" PER 72 CHARACTER LINE. 
Ee UY St Ohi ONOte 
YES 


THE REDUCED MODEL STATES ARE: 

1 2 3 ~ 5 6 1Z 13 14 15 43 heb 
46 47 48 49 50 aa 74 53 54 35 62 63 
65 66 67 68 69 70 Gal 72 a3 74 US 76 
78 iPS) 80 81 82 83 8y 85 86 87 88 89 
SA 2 93 94 $5 96 97, 98 


DO YOU, WESH TO CHANGE ANY OF THE SIGNIFICANT STATES? 


Ter ES LOR, NOM. 
NO 


ANALYSIS COMPLETE... YOUR REDUCED SYSTEM DATA HAS BEEN SAVED 
EN A. PLLE, NAMED “OFTMATR DATA". 


DO YOU DESIRE ANCTHER RUN? 


LPL. ES, SHE RESULTS OF THES RUN WILL 
REPLACE YOUR LAST CPTMATR DATA. 


Lee “YES OR “HOM. 
NO 


selesiea oo OPlnRey LS NOW TERMENATED. 5. cc c:s « 


feet —17.70/18.51 12:25:07 
RECORD OFF 
END RECORDING OF TERMINAL SESSICN 


The eigensystems for both systems were analyzed for 
Stability and other desireable characteristics. Bode piots 


were obtained for both input/output combinations and these 


Sy 


Fe 


45 
64 
aa 
90 


graphs for the full and reduced systems .are included as 
Figures 4.3 through 4.10. The input/output combination 
(1/1) shows the response of the system without the notch 
filter while the (1/2) input/output combination depicts the 
freguency response of the system after the addition of the 
notch,” £2iter: These figures clearly indicate excellent 
agreement between the full and reduced order systems in the 


frequency domain. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


See SS SS ——S— SS SS SS SS 


A. CONCLUSIONS 


The computational capabilities of OPTRED were tested by 
Bialyzing) two Significant large order systems. The, determi- 
Mation of the validity of both full order and reduced order 
results was greatly simplified by the eigensystem data for 
the X-29A provided by NASA-Edwards. As discussed by Laptas 
[Ref. 6], OPTSYSX failed to produce valid open-loop transfer 
function zeroes when analyzing very large systems. However, 
through modal analysis and with the aid of the NASA data, 
the extraneous zeroes were easily identified and the open 
loop data created by OPTSYSX was corrected. 

the results of) . OPTSYSX, and, ;OPTGRAPH.. indicate, a very 
strong Similarity between the full system and the reduced 
system generated by OPTRED both in stability and in the 
frequency domain. The combined analytic capabilities of 
OPTSYSX, OPTGRAPH and OPTRED provide a powerful tool for the 
controls designer for use in the analysis, development and 


the implementation of complex large order systems. 


B. RECOMMENDATIONS 


Based upon the results of this thesis the following 


areas are recommended for further study and research. 
1. Comparison of System Time Response 


Although the reduced model must closely resemble the 
full system in the frequency domain, it should also compare 
favorably in the time domain. OPTCALC, the time response 
program, requires several modifications to enable its use 


for analysis of systems as large as those presented in this 
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thesis. The addition of this frogram to the analysis tech- 
hique will provide a further means of verification of the 


accuracy and fidelity of the reduced order model. 


2. Interface with Other Analysis Programs 


— — —_—_  —— =—— —— 


The Naval Postgraduate SeEhoo's Aeronautical 
Engineering Department has developed a collection of anal- 
ysis programs including” OPITS¥sx, OPTGRAPH and OPTCALC. 
These programs exist ina single "module" which enabies the 
user to interactiveiy execute these programs during a single 
computer terminal session. The feasibility of adding OPTRED 
to this module should be investigated. If necessary, these 
programs would require re-dimensioning and could createa 
very large memory size requirement. Additions and modifica- 
tions would also be necessary in the OPTSYSX EXEC to enable 
the execution of OPTRED and then transfer execution to 
OPTSYSX for the analysis of the reduced model. This 
controlling EXEC must also manage the data files upon which 


the individual programs depend. 


3. Alternate Order Reduction ITechnigues 


Finally, other methods of order reduction should be 
investigated. Comparisons of reduced order models which are 
generated by different techniques would provide the controls 
designer with additional means ci reduced model verification 


and analysis. 
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